Effects of electron-phonon coupling range on the polaron formation 
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The polaron features due to electron-phonon interactions with different coupling ranges are in- 
vestigated by adopting a variational approach. The ground-state energy, the spectral weight, the 
average kinetic energy, the mean number of phonons, and the electron-lattice correlation function 
are discussed for the system with coupling to local and nearest neighbor lattice displacements com- 
paring the results with the long range case. For large values of the coupling with nearest neighbor 
sites, most physical quantities show a strong resemblance with those obtained for the long range 
electron-phonon interaction. Moreover, for intermediate values of interaction strength, the corre- 
lation function between electron and nearest neighbor lattice displacements is characterized by an 
upturn as function of the electron-phonon coupling constant. 
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I. INTRODUCTION 



In the last years effects due to strong electron-phonon (el — ph) interactions and polaronic signatures have been 
evidenced in several compounds, such as high-temperature cuprate superconductors,^: colossal magnetoresistance 
manganiteSf^i fuUerenes^^ carbon nanotubes^i and DNA^ This amount of experimental data has stimulated the study 
of different el—ph coupled systems. Among the most studied models there are the Holstein lattice model^ characterized 
by a very short-range (SR) el — ph interaction, and the Frohlich model^ that takes into account long-range (LR) 
el — ph couplings in polar compounds treating the dielectric as a continuum medium. Moreover, more realistic lattice 
interaction models including both SR and LR couplings have been recently introducedi^ 

Within the Holstein model a tight-binding electron locally couples to optical phonon modes. For intermediate el—ph 
couplings and electron and phonon energy scales not well separated, it has been found by numerical studies^iifliliiiSiii 
and variational approache3i^'iS*i& that the system undergoes a crossover from a weakly dressed electron to a massive 
localized polaronic quasiparticle, the small Holstein polaron (SHP). A variational approachi^*i& proposed by some of 
us and based on a linear superposition between the Bloch states characteristic of the weak and strong coupling regime 
is able to describe all the ground state properties with great accuracy. Moreover this method provides an immediate 
physical interpretation of the intermediate regime characterized by the polaron crossover. 

Recently a discrete version of the Frohlich model has been introduced in order to understand the role of LR coupling 
on the formation of the lattice polaron»ii Due to the LR interaction, the polaron is much lighter than the SHP with 
the same binding energy in the strong coupling regimefiS, Furthermore the lattice deformation induced by the electron 
is spread over many lattice sites giving rise to the formation of a large polaron (LP) also in the strong coupling 
regionjiSi Extending the variational approach previously proposed for the study of systems with local el — ph coupling, 
many properties have been studied by some of us as a function of the model parameters focusing on the adiabatic 
regime]^ Indeed there is a range of values of the el — ph coupling where the ground state is well described by a particle 
with a weakly renormalized mass but a spectral weight much smaller than unity. Furthermore, with increasing the 
strength of interaction in the same regime, the renormalized mass gradually increases, while the average kinetic energy 
is not strongly reduced. Finally due to the LR coupling a strong mixing between electronic and phononic degrees of 
freedom has been found even for small values of the el — ph coupling constant. 

Both the SR Holstein and the LR discrete Frohlich model can be described by a quite general Hamiltonian H 



H = -t c\cj -f Wo ^ a\ai -f aujQ ^ f{\Ri - Rj\)c\ci {a^ -f a]^ , (1) 

<i;j> i i,j 

where f{\Ri — Rj\) is the interacting force between an electron on the site i and an ion displacement on the site j. 
In Eq.(^ c| (ci) denotes the electron creation (annihilation) operator at site i, whose position vector is indicated by 
Ri, and the symbol <> denotes nearest neighbors (nn) linked through the transfer integral t. The operator al (ui) 
represents the creation (annihilation) operator for phonon on the site i, wq is the frequency of the optical local phonon 
modes, and a controls the strength of el — ph coupling. The units are such that h = 1. The Hamiltonian iQJ reduces 
to the Holstein model for 

f{\R.-R,\)^S^^j^^, (2) 
while in the LR case^^ the interaction force is given by 

_ 3 

f{\R,-R,\)^(\R,-R,\' + l) \ (3) 



if the distance \Ri — Rj\ is measured in units of lattice constant. 

In addition to the SR and LR case, in this work we analyze the properties of the system where the electron couples 
with local and nn lattice displacements. In this case the interaction force becomes 

fm-R,\)-ha. + '^ll^HMf (4) 

s 

where 5 indicates the nn sites. For all the couplings of Eqs. H2I3I4() the el — ph matrix element in the momentum 
space Mg- is 



^ L 



e^^-^™, (5) 
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with L number of lattice sites. Through the matrix element we can define the polaronic shift Ep 

and the coupling constant A = Ep/zt, with z lattice coordination number, that represents a natural measure of the 
strength of the el — ph coupling for any range of the interaction. Limiting the analysis to the one-dimensional case, 
the matrix element Mg- is reported in Fig. 1 for the SR, LR, and nn extended range {ER) couplings. In the SR case 
the coupling is constant as function of the transferred phononic momentum, while in the LR case the vertex is peaked 
around g ~ 0. With increasing the ratio ai/a, the ER interaction deviates from the constant behavior developing a 
peak around g ~ 0. Actually for the ratio ai/a — 0.3 the interaction vertex of the ER case is close to the behavior 
of the LR coupling. 

In this paper we adopt the variational approach previously proposedi^*iSi2£ for the study of systems with local 
and LR el — ph interactions in order to study the system with nn ER el — ph coupling. The aim is to investigate 
the crossover from SR to LR interactions in the ER model. The evolution of the ground-state spectral weight, the 
average kinetic energy, the mean number of phonons, and the electron-lattice correlation function with respect to the 
adiabaticity ratio loq /t and the el — ph coupling constant is discussed comparing the results with the local and LR 
case. For large values of the nn el — ph coupling many properties show a behavior similar to those obtained for the 
LR el — ph interaction. Regions of the model parameters are distinguished according to the values assumed by the 
spectral weight giving rise to a polaronic phase diagram. The transition line between the crossover and the strong 
coupling regime continuously evolves toward that of the LR case by increasing the coupling of the ER system. Finally, 
for intermediate values of this coupling, the correlation function between electron and nn lattice displacements shows 
an upturn with increasing the el — ph constant A. 

In section II the variational approach is reviewed, while in section III the results are discussed. 



II. VARIATIONAL WAVE FUNCTION 

In this section the variational approach is briefly summarized. Details can be found in previous worksi^ii^ 
The trial wave functions are translational invariant Bloch states obtained by taking a superposition of localized 
states centered on different lattice sites 

l€^>=;^E^''-'^"l4^(^")>, (7) 

-Rr. 

where 

l^fiR.) >= ■] 5:4)(i?„)4,jO > . (8) 

m 

In Eq. {7)) the apex i = w,s indicates the weak and strong coupling polaron wave function, respectively, |0 > denotes 
the electron and phonon vacuum state, and are variational parameters defining the spatial broadening of the 

electronic wave function. The phonon distribution functions h^^'^ (k) are chosen in order to reproduce polaron features 

in the two asymptotic limits^ 

In the intermediate regime the weak and strong coupling wave functions are not orthogonal and the off-diagonal 
matrix elements of the Hamiltonian are not zero. Therefore the ground state properties are determined by considering 
as trial state lip^ > a linear superposition of the weak and strong coupling wave functions 

/Al + Bl + 2ArBrSt 



where 
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and 




(11) 
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is the overlap factor of the two wave functions \tpj: > and >. In Eq.® Aj: and Bj: are two additional variational 

parameters which provide the relative weight of the weak and strong coupling solutions for any particular value of k. 
The variational minimization is performed extending the electron wave function up to fifth neighbors. 



III. RESULTS 



In this section we discuss ground state properties in the one-dimensional case for the different ranges of el — ph 
coupling. 

In Fig. 2(a) we report the polaron ground state energy as a function of the el — ph constant coupling A. The 
variational method recovers the perturbative results and improves significantly these asymptotic estimates in the 
intermediate region. In this regime the energy decreases with increasing the range of the e — ph coupling. Moreover, 
with increasing the range of the coupling, the crossover between the weak and strong coupling solution becomes less 
evident. Actually, as shown in Fig. 2(b), there are marked differences in the ratio B/A that is the weight of the strong 
coupling solution with respect to the weak coupling one. In the SR case the strong coupling solution provides all 
the contribution since the overlap with the weak coupling function is negligible. However, with increasing the range 
of the interaction, the weight of the weak coupling function increases and the polaronic crossover becomes smooth. 
Another quantity that gives insight about the properties of the electron state is the average kinetic energy K reported 
in Fig. 2(c) (in units of the bare electron energy). While in the SR case K is strongly reduced, in the LR case it is 
only weakly renormalized stressing that the self-trapping of the electron occurs for larger couplings with increasing 
the range of the interaction. Finally the mean number of phonons is plotted in Fig. 2(d). In the weak-coupling 
regime the interaction of the electron with displacements on different sites is able to excite more phonons. However, 
in the strong coupling regime there is an inversion in the roles played by SR and ER interaction. Indeed the SHP 
is strongly localized on the site allowing a larger number of local phonons to be excited. 

In addition to the quantities discussed in Fig. 2, other properties change remarkably with increasing the ratio ai/a. 
An interesting property is the ground state spectral weight Z, that measures the fraction of the bare electron state 
in the polaronic trial wave function. As plotted in Fig. 3(a), the increase of the el ~ ph coupling strength induces 
a decrease of the spectral weight that is more evident with increasing the range of the el — ph coupling. Only in 
the strong coupling regime the spectral weights calculated for different ranges assume similar small values. While for 
the local Holstein model Z — m/m*, as the ER case is considered, Z becomes progressively smaller than m/m* in 
analogy with the behavior due to the LR interaction. We have found that for the ratio ai/a = 0.3 there is a region 
of intermediate values of A where the ground state is described by a particle with a weakly renormalized mass but a 
spectral weight Z much smaller than unity. In Fig. 3(b) we propose a phase diagram based on the values assumed 
by the spectral weight making a comparison with that obtained in the SR and LR caseiiSiSS Analyzing the behavior 
of Z it is possible to distinguish different regimes, for example the crossover regime (0.1 < Z < 0.9) characterized by 
intermediate values of spectral weight and a mass not strongly enhanced for the ER case, and strong coupling regime 
{Z < 0.1) where the spectral weight is negligible and the mass is large but not enormous if the range of the coupling 
increases. With increasing the range of the interactions in the adiabatic case there is strong mixing of electronic and 
phononic degrees of freedom for values of A smaller than those characteristic of local Holstein interaction. Furthermore, 
entering the strong coupling regime, the charge carrier is still mobile and it does not undergo any abrupt localization. 
Only in the antiadiabatic regime the transition lines separating the crossover from the strong coupling regime tend 
to superimpose. 

Another important quantity associated to the polaron formation is the correlation function S{Ri) 



In Fig. 4(a) we report the correlation function S{Ri =0) at wo/t = 1 for several ranges of the el — ph interaction. In 
analogy with the behavior of the average number of phonons discussed in Fig. 2(d), the on-site correlation function 
is larger with increasing the range of the interaction in the weak coupling regime, but it becomes smaller as function 
of the coupling constant A in the strong coupling region indicating the the SHP is more effective in producing local 
lattice distortions. Actually in the ER and LR case, the lattice deformation is spread over nn or many lattice sites, 





(12) 
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respectively, giving rise to the formation of LP also in the strong coupling regime. Therefore it is interesting to analyze 
the behavior of the correlation function at nn sites. As reported in Fig. 4(b), in the SR case there is a minimum as 
function of A since the particle tends to localize on a single site with increasing the el — ph coupling. In the LR case 
the lattice distortion shows a decreasing behavior with increasing A indicating that the nearest neighbor contribution 
is always relevant. However for intermediate values of the ratio ai/a in the ER case, the correlation function shows 
an upturn as function of the coupling constant A. Actually, for small values of the coupling, this function tends to 
follow the behavior of the local interaction, but, with increasing the value of A, the coupling to the nn lattices is able 
to give deviations from the SR case. In fact the lattice deformation reaches a maximum, then begins to decrease 
following the behavior of the LR interaction. Therefore, as function of A, two different regimes in the correlation 
function can be evidenced. 



IV. DISCUSSION AND CONCLUSION 

In this paper we have extended a variational approach in order to study the polaronic ground-state features of a 
one dimensional el — ph model with coupling to local and nn lattice displacements. Many physical quantities such 
as the ground state energy and spectral weight, the average kinetic energy, the mean number of phonons, and the 
electron-lattice correlation function have been discussed making a comparison with the results obtained with SR and 
LR interactions. It has been possible to ascertain that most physical quantities are quantitatively equal to those 
obtained for the LR interaction as the el — ph coupling in the ER case is large. A polaronic phase diagram based on 
the values assumed by the spectral weight has been proposed. It has been shown that the transition lines between the 
crossover and the strong coupling regime continuously evolve toward that of the LR case by increasing the coupling 
of the ER system. The deviations of the ER case from LR case become evident only in quantities depending on 
distances larger than the lattice parameter, such as in the electron-lattice correlation function. At neighbor nearest 
sites for large values of the coupling, the ER interaction is able to reproduce the correlation function characteristic of 
the LR, while, at intermediate values of the ratio ai/a, the lattice deformation shows an upturn as function of the 
coupling constant A. 

Recently, a variational wave functionS^ has been proposed to study the polaron formation in Su-Schrieffer-Heeger 
(SSH) model where the electronic transfer integral depends on the relative displacement between nn sites. Unlike the 
original SSH model, the non-local electron-lattice coupling has been assumed to be due to the interaction with optical 
phonon modes. It has been shown that with this type of interaction the tendency towards localization is hindered 
from the pathological sign change of the effective next-nearest-neighbor hopping. Therefore it is not possible to reach 
the strong coupling regime where most properties obtained with the ER density-type el — ph coupling bear strong 
resemblance with those in the LR model. Only the coupling with acoustic phonons is able to provide a solution with 
localized behavior within the SSH modeli^ 

The variational approach for models with density-type el — ph coupling can be generalized to high dimensions, 
where it can still give a good description of ground state features^^*2i However, in order to reproduce with the ER 
interaction most physical quantities of the LR case, with increasing the dimensionality, it is important to include not 
only coupling terms at nn sites but also at next nearest neighbors. Actually it is necessary that the expansion of 
the coupling to near sites gives rise to an el — ph interaction vertex similar to that obtained in the LR case. Under 
these conditions the variational method is able to interpolate between the behavior of the SR case to the LR one 
with increasing the coupling of the interaction with close sites. 

Figure captions 

Fig.l The el — ph matrix element A/g (in units of awo/V^) for different ranges of the interaction as function of the 
momentum q (in units of tt). 

Fig. 2 The ground state energy i?o in units of cjo (a), the ratio B/A at k=0 (b), the average kinetic energy K in units 
of the bare one (c) and the average phonon number N (d) for t = as a function of the coupling constant 
A for different ranges of the el — ph interaction: SR (solid line), ER with ai/a = 0.05 (dash line), ER with 
aija = 0.1 (dot line), ER with aija — 0.2 (dash-dot line), ER with aijoL — 0.3 (dash-double dot line), LR 
(double dash-dot line). 

Fig. 3 (a) The ground state spectral weight at wo/t = 1 as a function of the coupling constant A for different ranges 
of interaction: SR (solid line), ER with aija ~ 0.05 (dash line), ER with ai/a = 0.1 (dot line), ER with 
aija — 0.2 (dash-dot line), ER with a\joL — 0.3 (dash-double dot line), LR (double dash-dot line). 
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(b) Polaron phase diagram for SR (solid line), ER with ai/a = 0.2 (dash-dot line) , ER with ai/a = 0.3 
(dash-double dot line), and LR (double dash-dot line) el — ph interaction. The transition lines correspond to 
model parameters such that the spectral weight Z = 0.1. 

Fig. 4 The electron-lattice correlation functions S{Ri = 0) (a) and S{Ri = S) (b) at wo/i = 1 for different ranges of 

the el — ph interaction: SR (solid line), ER with ai/a = 0.05 (dash line), ER with ai/a = 0.1 (dot line), ER 
with ai/a = 0.2 (dash-dot line), ER with ai/a = 0.3 (dash-double dot line), LR (double dash-dot line). 
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